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Title  of  Study:   Effects  of  P i nyon- Jun i per  Conversion  on  Watershed 

Values  in  Utah 
Object  i ves : 

A.  To  determine  the  water  budget  of  natural  stands  of  pinyon- 
juniper  and  adjacent  areas  which  have  been  cleared  and/or  seeded. 

B.  To  determine  the  effects  of  vegetation  conversion  on  soil 
physical  properties  and  soil  stability. 

C.  To  ecologically  evaluate  sites  before  and  after  as  to 
composition  and  production  of  vegetation. 

D.  To  evaluate  the  economics  of  conversion  practices  in  terms 
of  the  watershed  values  and  multiple  use  relations. 

E.  To  obtain  data  necessary  for  determination  of  hydrologic 
soil  cover  complexes  on  the  study  sites. 

Introductory  Comment:  This  report  is  concerned  with  additional  data 
analysis  and  compilation  which  has  resulted  since  the  project  report 
dated  April  1,  1971.  As  before,  the  report  will  provide  information 
to  supplement  previous  reports  as  well  as  indicate  progress  to  date. 

Runoff  Plot  Studies:   No  runoff  events  occurred  at  either  site 
during  1971 • 

Tables  1,  2  and  3  give  precipitation  data  from  the  Blanding  study 
site.   These  tables  can  be  summarized  on  one  word  dry. 

Tables  k,    5,  and  6  give  rainfall  data  from  the  Milford  study 
site.   There  was  a  complete  absence  of  runoff-producing  storms. 


Table    l  .       Precipitation  data  from  8-inch  recording  gages  at  Blanding 
study  site,  1971. 


Date 


Windrow  Area 


Total  Rainfall   (inches) 


Hebris-in-Place 


6-11-71 

6-12-71 

7-21-71 

8-5-71 

8-7-71 

8-17-71 

8-19-71 

8-20-71 

8-21-71 

8-23-71 

8-24-71 

8-25-71 

8-26-71 

8-28-71 

8-30-71 

9-30-71 

10-1-71 

10-11-71 


Start 

0.17 
0.26 
0.03 
0.20 
0.25 
0.51 
0.40 
0.03 
0.04 
0.45 
0.45 
0.18 
0.18 
no  record 
no  record 


Start 

0.21 

0.24 

0.02 

0.20 

0.33 

0.53 

0.24 

0.04 

0.05 

0.40 

0.45 

0.23 

0.21 

0.55 

0.38 


OFF 
Storage  gage  charged 


Table    2.     Precipitation  data  from  8-inch  nonrecording  gages  at  Blanding, 
debris-in-place  area,  1971. 


Total  Rainfall (inches) 


Date  Gage  A  Gage  B 


6-11-71   to  6-28-71  0.00  0.00 

6-28-71   to  7-10-71  0.00  0.00 

7-10-71   to  7-24-71  0.22  0.23 

7-24-71   to  8-7-71  0.24  0.23 

8-7-71   to  8-20-71  0.53  0.43 

8-20-71   to  9-3-71  2.23  2.16 

9-3-71   to  9-14-71  0.00  0.00 

9-14-71   to  10-11-71  0.93  0.90 


Table    3    .Precipitation  data  from  8-inch  nonrecording  gages  at  Blanding, 
windrow  area,  1971 . 

Total  Rainfall (inches) 

Date  Gage  A  Gage  B 


10-26-70  to  4-3-71  2.20* 

4-4-71   to  6-11-71  0.90* 

6-12-71   to  6-28-71  0.00  0.00 

6-28-71   to  7-9-71  0.00  0.00 

7-9-71   to  7-24-71  0.18  0.17 

7-24-71   to  8-7-71  0.29  0.25 

8-7-71   to  8-20-71  0.48  0.45 

8-20-71  to  9-3-71  2.30  2.05 

9-3-71   to  9-14-71  0.00  0.00 

9-14-71   to  10-11-71  0.83  0.90 


*Single  storage  gage  operated  during  this  period. 


Table      h.       Precipitation  data  from  8-inch  recording  gages  at  Afford 
study  site,  1971. 

Total  Rainfall   (inches) 
Date Windrow  Area Debris-in-Place 

6-8-71  Start 

6-9-71 

6-11-71 

7-15-71 

7-16-71 

7-17-71 

7-19-71 

7-20-71 

7-21-71 

7-22-71 

7-30-71 

8-3-71 

8-5-71 

8-6-71 

8-8-71 

8-14-71 

8-15-71 

8-16-71 

8-17-71 

8-18-71 

8-20-71 

8-21-71 

8-26-71 

8-27-71 

8-28-71 

9-30-71 

10-1-71 

10-6-71 


0.14 

0.19 

0.08 

0.08 

0.05 

0.04 

0.13 

0.15 

0.09 

0.08 

0.05 

0.04 

0.12 

0.12 

0.25 

0.22 

0.07 

0.06 

0.27 

0.28 

0.15 

0.14 

0.06 

0.07 

0.03 

0.03 

no  record 

0.04 

0.10 

0.09 

0.06 

0.03 

0.18 

0.13 

0.70 

0.74 

0.41 

0.46 

0.10 

0.08 

0.22 

0.26 

0.25 

0.29 

0.15 

0.18 

0.25 

0.18 

0.45 

0.45 

OFF 

Storage 

gage  charged 

Table    5.     Precipitation  data  from  8-inch  nonrecording  gages  at  Mi  1  ford , 
debris-in-place  area,  1971. 

Total  Rainfall (inches) 
Date 


6-9-71  to  6-22-71 
6-22-71  to  7-6-71 
7-6-71  to  7-20-71 
7-20-71  to  8-3-71 
8-3-71  to  8-17-71 
8-17-71  to  9-6-71 
9-6-71  to  9-16-71 
9-16-71   to  10-6-71 


iage  A 

Gage  B 

Gage  C 

0.15 

0.10 

0.22 

0.00 

0.00 

0.00 

0.37 

0.32 

0.36 

0.48 

0.45 

0.45 

0.79 

0.77 

0.85 

2.05 

2.00 

2.35 

0.00 

0.00 

0.00 

0.78 

0.62 

Table    6.     Precipitation  data  from  8-inch  nonrecording  gages  at  Mil  ford, 
windrowed  area,  1971. 


Total  Rainfall   (inches) 


Date  Gage  A  Gage  B  Gage  C 

10-24-70  to  4-10-71 
4-11-71  to  6-8-71 
6-8-71  to  6-22-71 
6-22-71  to  7-6-71 
7-6-71  to  7-20-71 
7-20-71  to  8-3-71 
8-3-71  to  8-17-71 
8-17-71  to  9-6-71 
9-6-71  to  9-16-71 
9-16-71   to  10-6-71 


1.10* 

3.00* 

0.15 

0.03 

0.10 

0.00 

0.00 

0.00 

0.35 

0.31 

0.31 

0.53 

0.56 

0.54 

0.87 

1.00 

1.05 

1.95 

1.95 

1.80 

0.00 

0.00 

0.00 

0.89 

0.84 

0.82 

*Single  storage  gage  operated  during  this  period. 


Influence  of  Graz  i  ng  on  Se lected  Hydrol og  ic  Parameters 
of  P  inyon- Jun  iper  Range  land 

Grazing  is  the  major  use  of  the  pinyon  (P  inus  spp.)  and  juniper 
( Jun  i  perus  spp.)  vegetation  type  in  Utah,  Arizona,  Colorado  and  New  Mexico 
(Aro  1971).   Hydrologic  consequences  of  grazing  are  determined  by  the 
kind  and  degree  of  modifications  made  in  the  vegetal  cover,  surface 
and  subsurface  soil  conditions,  and  intensity  of  rainfall  characteristics 
of  the  area  in  which  the  changes  occurred  (Rauzi  and  Kuhlman  19&1). 

Woodbury  (19^7),  Mason  (1963),  and  Aro  (1971)  reported  that 
overgrazing  of  pi nyon- juniper  rangeland  has  resulted  in  reduced  forage 
production  and  increased  extent  and  density  of  the  tree  species.   Managers 
of  pinyon- jun iper  range  seek  to  convert  these  woodlands  into  more  productive 
grasslands.   Public  agencies  supervising  these  vegetation  manipulation 
programs  are  also  interested  in  protection  or  improvement  of  wildlife, 
aesthetic,  and  watershed  values  of  the  sites.   This  study  involves  the 
continuing  effort  to  evaluate  hydrologic  effects  of  converting  pinyon- 
juniper  woodland  to  open  grassland. 


OBJECTIVES: 

Objectives  of  this  study  are  (1)  to  determine  the  effects  of  grazing, 
when  compared  to  areas  on  which  grazing  has  been  excluded,  on  selected 
hydrologic  parameters  of  chained,  debr is- in-place ;  chained,  debris- 
windrowed;  and  unchained  pinyon- jun i per  sites  and  (2)  utilize  these 
measurements i n  developing  guidelines  for  grazing  management  of  pinyon- 
juniper  rangeland  that  will  protect  and/or  improve  the  hydrologic  value 
of  the  watershed.   Of  particular  interest  are  the  following  items: 
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1.  Influence  of  grazing  on  infiltration  and  sediment  production 
rates,  especially  as  related  to  surface  soil  changes  and 
vegetal  modifications, 

2.  Changes  in  infiltration  and  sediment  production  rates 
as  a  function  of  time  since  grazing  has  been  excluded, 

3.  Develop  multiple  regression  models  for  predicting  infiltration 
and  erosion  rates  of  chained  and  unchained  pi nyon- jun i per 
rangeland  subjected  to  various  grazing  situations. 

The  following  treatments  are  being  investigated: 

A.  Chained,  debr is- in-place . 

1.  Grazing  not  excluded 

2.  Grazing  excluded  1 967 

3.  Grazing  excluded  1969 

4.  Grazing  excluded  1971 

B.  Chained,  debr is-wi ndrowed 

5.  Grazing  not  excluded 

6.  Grazing  excluded  1967 

7.  Grazing  excluded  1971 

C.  Not  chained 

8.  Grazing  not  excluded 

9.  Grazing  excluded  1967 

10.  Grazing  excluded  I969 

11.  Grazing  excluded  1971 

RESEARCH  PLANS: 

Initial  hypothesis  of  this  study  is  that  grazing  of  chained  and 
unchained  pi nyon-jun iper  rangeland  will  dffect  hydrologic,  vegetative,  and 
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edaphic  parameters  of  the  sites  and  these  effects  will  vary  with  the 
length  of  time  grazing  has  been  excluded  and  with  the  particular 
treatment  imposed.   Expected  effects  attributable  to  grazing  are 
removal  of  above  ground  plant  material;  and  therefore,  decreased  soil 
protection;  compaction  of  the  soil  surface  by  trampling;  altered 
infiltration  capacities;  and  accelerated  erosion.   Elimination  of  grazing 
is  expected  to  eliminate  or  reverse  these  hydrolog ical ly  detrimental 
effects.   The  presently  unknown  time  sequence  of  recovery  is  also 
being  investigated. 

Study  area.   Previously  chained  and  seeded  sites  near  Blanding, 
Utah  were  chosen  for  study.   The  area  is  approximately  in  the  center 
of  the  pinyon-jun i per  distribution  range  and  the  availability  of  closely 
adjacent  chained,  debr is- in-place  ;  chained,  debr is-wi ndrowed ;  and  unchained 
sites  with  grazing  exclosures  previously  established  for  one,  two,  and 
four  years  made  this  a  highly  desirable  study  area. 

The  areas  surrounding  the  grazing  exclosures  are  grazed  from  May  1 
to  June  15  and  October  1  to  November  1  each  year  with  cattle.   Management 
of  the  herd  does  not  allow  estimation  of  the  grazing  intensity. 

Basic  research  procedure.   Runoff  and  erosion  was  artifically 
induced  from  small  plots  by  simulating  high  intensity  (3  in/hr  or  greater) 
rainfall  with  the  Rocky  Mountain  Inf i 1 trometer  (Dortignac  1951).   The 
infiltration  rate  (rainfall  appl ied- runoff) ,  sediment  production,  and 
selected  vegetation  and  edaphic  parameters  were  measured  from  each 
plot  sampled.   Data  were  collected  immediately  after  the  spring  grazing 
period  and  after  two  months  of  rest  from  grazing  pressure  during  1971. 
Approximately  231  i nf i 1 trometer  plots  were  sampled  during  each  period. 
Field  procedures  were  similar  to  those  described  by  Williams  (196'i). 
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One-third  of  the  plots  (not  grazing  exclosures)  were  clipped  to 
simulate  grazing  occurring  on  adjacent  grazed  areas.   Another  one-third 
of  the  ungrazed  plots  were  completely  denuded  to  approximate  the  bare 
soil  condition  encountered  in  the  unchained  areas.   These  plots  provide 
information  on  hydrologic  effects  associated  with  forage  removal  but 
with  trampling  eliminated.   The  remaining  plots  received  no  pre- inf i 1 trat ion 
test  clipping  treatment.   Vegetation  and  soil  surface  cover  characteristics 
were  determined  for  all  plots  after  these  clipping  treatments  were 
imposed . 

Following  the  infiltration  and  erosion  measurements,  vegetative 
material  present  within  each  inf i 1 t raometer  plot  was  clipped  at  ground 
level,  oven  dried,  and  weights  converted  to  pound/acre  of  forage. 

One-and  3— inch   deep  (3  inches  diameter)  soil  cores  were  collected 
from  within  the  plots  and  hydraulic  conductivity,  bulk  density,  and 
porosity  determined  for  each  core  (Hoover,  Olson,  and  Metz  195M- 
These  parameters  were  determined  for  each  sampling  period.   In  addition, 
following  the  second  sampling  period  soil  texture  and  percent  water 
stable  aggregates  were  determined  for  the  two  soil  cores  (Bouyoucos 
1929  and  1962) . 

Stepwise  multiple  regression  modeling  techniques  will  be  utilized 
to  explain  variation  in  infiltration  and  erosion  rates  (dependent 
or  'y'  variables)  using  the  hydrologic, edaph ic ,  and  vegetative  parameters 
as  independent  or  'x1  variables. 

RESEARCH  ACCOMPLISHMENTS: 


All  information  described  above  was  collected  during  the  1971 
sampling  periods;  however,  data  analysis  is  not  complete. 


/2 


Maximum  infiltration  rates  for  each  treatment  and  sampling  period 
are  indicated  in  Figures  /  to  (p  .       In  general,  infiltration  rates  were 
lowest  on  grazed  plots,  with  increasing  rates  exhibited  as  the  period 
of  protection  from  grazing  increased.   Grazed  plots  recorded  lower 
water  intake  rates  than  the  plots  receiving  the  pre-run  clipping 
treatment.   Therefore,  reduced  infiltration  rates  found  on  grazed  areas 
may  not  be  attributable  to  forage  removal,  but  rather  to  soil  compactic 
by  the  an  ima 1  . 

The  June-July  sampling  period  (Figures  /  ,  £  ,   3  )    exhibited  a 
higher  minimum  infiltration  rate  than  the  August-September  period 
(Figures  $-   ,  S"  ,  £     )  .      This  difference  is  probably  due  to  a  seasonal 
change  in  infiltration  capacity.   This  late  season  reduction  in  water 
intake  rates  is  of  considerable  practical  importance  because  high 
intensity  thunderstorms  commonly  occur  during  this  period. 

FUTURE  RESEARCH  PLANS: 


An  artifical  trampling  treatment  will  be  included  in  the  study 
during  the  June,  1972,  sampling  period.   This  treatment  will  consist 
of  subjecting  randomally  selected  points  (covering  approximately  one 
half  of  the  plot  area)  to  a  compacting  force  of  2h   pounds  per  square 
inch.   This  force  is  the  static  load  pressure  calculated  for  cattle 
by  Lull  (1959).   This  compaction  technique  will  be  applied  to  one-third 
of  the  ungrazed  plots  which  have  also  been  clipped  to  a  simulated 
grazing  intensity.   This  treatment  should  approximate  conditions  occurring 
on  grazed  plots.   Another  one-third  of  the  ungrazed  plots  will  be  compacted 
but  not  clipped.   These  plots  should  provide  data  on  hydrologic  effects 
associated  with  trampling  but  with  forage  removal  eliminated.   The 
remaining  plots  will  not  be  disturbed  prior  to  the  infiltration  test. 
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Soil  Moisture  Studies: 

Milford.   Soil  moisture  data  are  now  available  for  two  years  from  the 
Milford  site.   Figures     _(windrow  vs.  control)  and  $   (debris-left-in- 
place  vs.  control)  represent  average  inches  of  water  per  foot  of  soil 
profile  on  various  sampling  dates.   The  windrowed  treatment  had  consistantly 
greater  soil  moisture  values  per  foot  of  soil  profile  than  did  respective 
controls.   The  differences  in  water  per  foot  of  soil  profile  between  the 
debris-in-place  treatment  and  respective  controls  were  uniformily  greater 
than  the  windrowed  treatment  and  its  respective  control. 

Regardless  of  treatment,  the  greatest  differences  in  water  per  foot 
of  soil  profile  occurred  during  March  through  April.   The  greatest 
accumulation  of  soil  moisture  was  in  early  spring. 

Figures   /  (windrow  vs.  control)  and  ft      (debris-in-place  vs.  control) 
represent  average  inches  of  water  per  foot  of  soil  profile  averaged  over  two 
years.   The  average  amount  of  water  is  least  in  the  first  foot  of  soil  profile. 
Average  soil  moisture  at  the  2-,  3-,  and  U-foot  depths  is  relatively  uniform. 
The  greatest  difference  in  soil  moisture  between  the  chaining  treatments 
and  their  controls  is  in  the  last  two  feet  of  soil  profile. 

Blanding.   Soil  moisture  data  are  available  for  one  year  from  the 
Blanding  site.   Figure  J  I     represents  average  inches  of  water  per  foot  of 
soil  profile  for  the  three  treatments  on  various  sampling  dates.   The 
debris-in-place  treatment  had  the  most  water  per  foot  of  soil  profile  over 
the  whole  year.   There  is  a  slight  peak  in  April,  which  probably  indicates 
the  time  of  maximum  accumulation  of  soil  moisture. 

Figure  IvL  represents  average  inches  of  water  per  foot  of  soil  profile 
for  each  treatment,  averaged  over  one  year.   As  at  Milford,  the  least  averatir<' 
amount  of  soil  moisture  is  in  the  first  foot  of  soil  profile.   From  2  to  5 
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feet  the  soil  moisture  is  relatively  uniform.  The  greatest  difference  in 
soil  moisture  among  the  different  treatments  is  in  the  last  three  feet  of 
soil  profile. 
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Sap  Velocity  Studies 

A  study  on  sap  flow  patterns  of  Utah  juniper  (Juniperus  osteosperma) 
and  pinyon  pine  (Pinus  edulis)  was  done  near  Blanding,  Utah,  over  a  two 
year  period.   The  object  of  the  study  was  to  compare  sap  velocities  of  the 
two  species.   Soil  moisture  and  other  environmental  factors  (wind, 
temperature,  relative  humidity  and  incoming  short  wave  radiation)  were 
taken  during  the  study.   Sap  flow  was  taken  three  to  four  days  per  month 
during  the  summer  and  one  day  or  less  per  month  during  the  rest  of  the 
year. 

The  data  were  statistically  analyzed  and  results  indicated  that 
soil  moisture  availability  was  the  major  factor  in  controlling  which 
species  had  the  highest  sap  velocities.   During  1970,  the  Blanding  area 
had  an  above  average  year  in  precipitation.   Pinyon  pine  trees  had  a 
significantly  higher  sap  velocity  than  the  Utah  juniper  trees  from  August 
until  the  end  of  November.   From  December  1970  until  September  1971,  the 
precipitation  was  way  below  normal  for  the  same  area.   There  was  no  significant 
difference  in  sap  velocities  between  the  two  species  of  trees  from  December 
until  May.   When  the  water  deficit  was  the  greatest  from  June  to  July, 
Utah  juniper  trees  had  a  significantly  higher  sap  velocity  than  the 
pinyon  pine  trees. 

It  is  concluded  that  when  soil  moisture  is  readily  available  the 
pinyon  pine  trees  has  the  significantly  highest  sap  velocity  of  the  two 
species.   Utah  juniper  trees  have  significantly  higher  sap  velocities  than 
pinyon  pine  trees  when  soil  moisture  is  not  readily  available.   A  possible 
explanation  for  the  different  responses  of  the  two  species  is  the  different 
physiological  characteristics  of  the  leaves  and  different  rooting 
characteristics . 


'"I 


Ran  velocity  was  no^  strongly  related  to  green  "biomass  of  either  tree 
species.   Figure  /3   relates  green  biomass  (both  species)  to  diameter  at 
the  one  foot  height  (r  =  .95*0  and  Figure  \*v   shows  the  relation  of  sap  velocity 
(for  both  species)  to  tree  diameter  at  the  one  foot  level.   This 
indicates  that  sap  velocities  are  apparently  independent  of  tree  diameter  an<.\ 
green  biomass,  meaning  that  total  water  use  is  a  function  of  conducting 
area  of  the  tree  stem,  as  well  as  sap  velocity. 
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Vegetat  ion  Stud  ies ; 

Tables  8  and  9  give  tree,  shrub,  and  ground  cover  on  debr i s- i n-p lace 
and  windrowed  runoff  plots,  respectively,  for  the  Milford  site,  1971. 

Tables  10  and  11  give  the  same  information  for  debr i s- i n-p lace 
and  windrowed  runoff  plots,  respectively,  for  the  Blanding  site,  1971- 


33 


TABLE  $       Tree,  shrub,  and  ground  cover  (percent)  on  debris-in-place  runoff 
plots  at  Milford  study  site,  September,  1971. 


Plot 


Transect   

Number  1/  Trees 


Percent  Cover 


Shrub 


Ground  2/ 


Debris-in- 
place  #1 


19  ft 


33  ft 


55  ft 


7k   ft 


(Mean) 


Artr 

9 

BG  2 

Chvi 

2 

L  30 

Agcr   6 

EP  58 

Misc   1 

Sphaeralcea  spp  3 

Artr 

1+ 

BG   0 

Chvi 

10 

L  22 

Agcr  8 

EP  68 

Chvi  10 

Sphaeralcea  spp   2 

Artr 

11 

BG  0 

Chvi 

1 

L  3^ 

Agcr  8 

EP  57 

Sphaeralcea  spp  1 

Artr 

3 

BG  1 

Chvi 

2 

L  29 

Agcr  8 

EP  58 

Sphaeralcea  spp  h 

Artr 

7 

BG  1 

Chvi 

h 

L  28 

Agcr  8 

EP  6l 

Sphaeralcea  spp  2 

Debris-in- 
place  #2 


19  ft 


33  ft 


Artr  9 

BG   3 

Chvi  11 

L  39 

Agcr   1 

EP  55 

Phlox  spp  2 

Artr   3 

BG   3 

Chvi  13 

L  38 

Agcr  7 

EP  U8 

Orhy  2 

Phlox  spp  1 

Sphaeralcea  spp  1 


Table  %     (contd) 
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■lot 


Transect        

Number     1/      Trees 


Percent   Cover 


Shrub 


Ground  2/ 


55   ft 


7k  ft 


(Mean)    x 


Chvi 

in 

BG   2 

L  37 

Agcr  ] 5 

EP  III 

Penstemon  spp  2 

Eriogonum  spp  2 

Sphaeralcea  spp   1 

Artr 

7 

BG  2 

Chvi 

8 

L  36 

Agcr  12 

EP  U9 

Lupine  spp  1 

Artr 

5 

BG  3 

Chvi 

11 

L  38 

Agcr   9 

EP  1+8 

Phlox  spp   1 

Lupine  spp   1 

Debris-in- 
place   #3 


19  ft 


33  ft 


55  ft 


lh  ft 


Juos   5 


Artr 

2 

BG 

k 

Chvi 

2 

L 

6o 

Agcr 

IT 

EP 

15 

Lupine  spp 

3 

Eriogonum  spp 

1 

Artr 

6 

BG 

12 

Chvi 

2 

L 

U3 

Agcr 

27 

EP 

13 

Eriogonum  spp 

1 

Penstemon  spp 

1+ 

Artr 

5 

BG 

3 

Chvi 

5 

L 

31 

Agcr 

18 

EP 

k6 

Eriogonum  spp 

1 

Sphaeralcea  spp 

1 

Artr 

1 

BG 

8 

Chvi 

8 

L 

39 

Agcr 

1  3 

EP 

37 

Penstemon  spp 

1 

Eriogonum  spp 

1 

Sphaeralcea  spp 

1 

,.?<r 


Table  $    (contd) 


Transect 

Perc 

ent  Cover 

Number  1/  Trees 

Shrub 

Ground 

2/ 

(Mean)   5c      Juos   1 

Artr  h 

BG 

6 

Chvi  k 

L 

Agcr 

EP 

k2 
19 
29 

Eriogonum  spp 

Lupine  spp 

Sphaeralcea  spp 

Penstemon  spp 

1 
1 
1 
1 

Depris-in- 

place  ffk              19  ft     Juos  15 

Artr  2 

BG 

1 

Chvi  10 

L 

Agcr 

EP 

Phlox  spp 

Sphaeralcea  spp 

hi 

18 

31 

2 

1 

33  ft         0 

Artr  lk 

BG 

1 

Chvi   7 

L 

Agcr 

EP 

Phlox  spp 

Sphaeralcea  spp 

31 
18 
kk 

1 
5 

55  ft         0 

Artr  1 

BG 

0 

Chvi   1 

L 

Agcr 

EP 

9 
15 
61 

Orhy 

1 

. 

Sphaeralcea  spp 

8 

Ik   ft          0 

Artr  3 

BG 

L 

Agcr 

EP 

Sphaeralcea  spp 

0 
10 

16 

6U 

8 

Eriogonum  spp  1 

Lupine  spp 

1 

(Mean)  x       Juos  k 

Artr   5 

BG 

1 

Chvi  k 

L 

Agcr 

EP 

25 
18 
51 

Phlox  spp  1 
Sphaeralcea  spp  k 


Table  %   (contd) 


3C, 


Plot 


Transect    

Number  1/  Trees 


Percent  Cover 


Shrub 


Ground  2/ 


Debris-in- 
place  #5 


19  ft 


(Mean)  x 


Pied  8 


33  ft     Pied  h 


55  ft     Pied  3 


7^  ft     Juos  1 


Pied  h 


Artr 

5 

BG  0 

Chvi 

10 

L  3h 

Agcr  12 

EP  51* 

Artr 

18 

BG   1 

Chvi 

6 

L  33 

Agcr  13 

EP  51 

Lupine  spp   1 

Eriogonum  spp  1 

Artr 

16 

BG  0 

Chvi 

10 

L  20 

Agcr  6 

EP  73 

Sphaeralcea  spp  1 

Artr 

k 

BG  0 

Chvi 

Ik 

L  20 

Agcr   8 

EP  71 

Eriogonum  spp  1 

Artr 

11 

BG  0 

Chvi 

10 

L  27 

Agcr  10 

EP  62 

Misc   1 

1/  Line  transects  across  runoff  plots  at  indicated  distances 
measured  from  top  of  plot. 

2/   BG  =  Bare  ground 
L  =  Litter 
Agcr  =  Agropyron  cristatum 

EP  =  Erosion  Pavement 
Pied  =  Pinus  edulis 
Juos  =  Juniperus  osteosperma 
Chvi  =  Chrysothamnus  viscidiflorus 
Orhy  =  Oryzopsis  hymenoides 


Table  /.   Tree,  shrub,  and  ground,  cover  (percent)  on  windrowed  runol 
plots  at  the  Mil ford  study  site,  September,  1971. 


3? 


Plot 


Transect 
Number  1/ 


Percent  Cover 


Trees 


Shrub 


Ground— 


2f 


Windrow  #1    19  ft 


33  ft 


55  ft. 


'fk   ft. 


(Mean)   x 


Artr  1 


BG 
L 
Agcr 
EP 

Rock 
Penst 


0 
I* 

30 

6k 

1 

emonl 


Artr  2 

BG 

3 

L 

6 

Agcr 

26 

EP 

6o 

Rock 

U 

Misc  . 

l 

0 

BG 

0 

L 

k 

Agcr 

29 

EP 

6k 

Rock 

D 

Artr  1 

BG 

6 

L 

10 

Agcr 

29 

EP 

53 

Rock 

2 

BG 
L 
Agcr 
EP 
Rock 
Penst 


2 

6 

29 

6o 

2 
emonl 


Windrow  #2    19  ft 


33  ft 


55  ft 
7k   ft 


(Mean)  x 


BG  70 

L  k 

Agcr  2k 

Lupine 


BG 

6l 

L 

16 

Agcr 

23 

BG 

6^5 

L 

9 

Agcr 

26 

BG  57 

L  12 

Agcr  30 

Misc.  1 

"bg  oT 

L  10 

Agcr  26 

Lupine  1 


Table  /  .   Continued, 


36 


Plot 


Transect 
Number—' 


Percent  Cover 


Trees 


Shrubs 


Groundi 


Windrow  #3    19  ft. 


33  ft 


55  ft 


lh   ft, 


(Mean)   x 


BG 
L 
Agcr 
EP 

Penstemon 
Misc  . 


BG 
L 
Agcr 
EP 
Rock 


BG 
L 
Agcr 
EP 
Rock 


BG 
L 
Agcr 
EP 
Penstemon 


BG 
L 
Agcr 
EP 

Rock 
Penst 


3 

2 

23 

TO 
1 
1 


1 

0 

21 

69 
i 

"T 

7 
23 
65 

1 


18 
18 
23 
36 
5 


7 

9 

22 

60 

1 

emon   1 


Windrow  ffh  19  ft. 


33  ft 


55  ft. 


7^  ft 


Artr  1 


Artr  1 


BG 
L 
Agcr 
EP 

Lupine 
Penstemon 


BG 
L 
Agcr 
EP 

Lupine 
Penstemon 


BG 
L 
Agcr 
EP 
Lupine 


BG 

L 
Agcr 
EP 
Lupine 


20 

1 

25 

hk 

9 

1 


36 
10 
21 

31 
1 

1 


32 
3 
2k 
31 
10 


57 

12 

23 

0 

8 


Table  /  .   Continued. 


39 


Plot 


Transect 


Numberi-' 


1/ 


Trees 


Percent  Cover 


Shrubs 


Ground— 


2t 


(mean)  x~ 


Artr  1 


BG 
L 
Agcr 
EP 

Lupine 
Pensteraon 


37 

6 

23 

26 

7 

1 


Windrow  #5    19  ft. 


33  ft. 


55  ft, 


Ik   ft, 


(Me 


an   x 


BG 
L 
Agcr 
Lupine 


BG 
L 
Agcr 


66 
1 

31 
2 


71 

2 

27 


BG 

55 

L 

ik 

Agcr 

30 

Rock 

1 

BG 

70 

L 

5 

Agcr 

22 

Rock 

3 

BG 

65 

L 

6 

Agcr 

27 

Lupine 

l 

Rock 

1 

—Line  transects  across  runoff  plots  at  indicated  distances  measured  from 
top  of  plot. 

2/BG  =  Bare  Ground 

L  =  Litter 
Agcr  =  Agropyron  cristatum 
EP  =  Erosion  Pavement 


4i 


TABLE  //'/Tree,  shrub,  and  ground  cover  (percent)  on  debris-in-place  runoff 
plots  at  Blanding  study  site,  September,  1971 • 

p         Transect    Percent  Cover 

Number  1/   Trees     Shrub  Ground  2/ 

Debris-in- 
place  #1      19  ft     Juos  18        0 


33  ft 


55  ft 


Ik   ft     Juos  22 


(Mean)  x       Juos  10 


Debris-in- 
place  #2      19  ft     Juos   U 


33  ft 


55  ft 


7^  ft     Juos  9 


BG 

10 

L 

83 

Agcr 

6 

Saka 

1 

BG 

32 

L 

35 

Agcr 

27 

Saka 

5 

Ast 

er  spp 

1 

BG 

33 

L 

1+7 

Agcr 

20 

BG 

51 

L 

1+0 

Agcr 

7 

Saka 

1 

Ast' 

er  spp 

1 

BG 

32 

L 

51 

Agcr 

15 

Ast> 

sr  spp 

1 

Saka 

1 

BG 

5^4 

L 

19 

Agcr 

27 

BG 

1+2 

L 

50 

Agcr 

8 

BG 

36 

L 

51 

Agcr 

13 

BG 

3h 

L 

50 

Agcr 

12 

Saka 

k 

Table  /O    (contd) 


4-/ 


Plot 


Transect 

Number  1/   Trees 


Percent  Cover 


Shrub 


Ground  2/ 


(Mean)  x 


Juos   3 


BG  111 

L  k2 

Agcr  16 

Saka  1 


Debris-in- 

place  ~'3 


19  ft 


33  ft 


55  ft 


(Mean)  x 


Juos 


lk   ft     Juos  17 


Juos    6 


BG 

kk 

L 

32 

A 

Agcr 

2k 

BG 

Ik 

L 

lk 

Agcr 

12 

BG 

30 

L 

53 

Agcr 

17 

BG 

k5 

L 

33 

Agcr 

22 

Sphaeral 

cea  spp 

2 

BG 

33 

L 

kQ 

Agcr 

19 

Debris-in- 
place   ttk 


19   ft 


33  ft 


55  ft 


Ik   ft 


BG 

k2 

L 

k3 

Agcr 

13 

Saka 

2 

BG 

30 

L 

3k 

Agcr 

30 

Saka 

2 

Sphaeralc 

ea  spp 

2 

BG 

52 

L 

27 

Saka  2 

Agcr 

18 

Sphaeralcea  spp 

1 

BG 

)i 

L 

88 

Agcr 

5 

Sphai  ira  I  cea  spp 

3 

Table  /O     (contd) 


&   1 


Plot 


Transect 


Trees 


Percent  Cover 


Shrub 


Ground  2/ 


(Mean)  x 


BG  32 

L  1+8 

Agcr  l6 

Saka  2 

Sphaeralcea  spp  2 


Debris-in- 
place  #5      19  ft 


33  ft 


55  ft     Juos  6 


Ik   ft 


(Mean)  x 


Juos  1 


BG 

51 

L 

36 

Agcr 

13 

BG 

Uo 

L 

kh 

Agcr 

16 

BG 

29 

L 

50 

Agcr 

5 

Ast< 

sr  spp 

2 

Saka 

6 

Opunt: 

La  spp 

2 

Sphaeralcea  spp 

6 

BG 

61 

L 

22 

Agcr 

IT 

BG 

^5 

L 

1+0 

Agcr 

13 

Saka 

1 

Sphaeralcea  spp 

1 

1/  Line  transects  across  runoff  plots  at  indicated  distances  measured 
from  top  of  plot. 


2/    BG  =  Bare  ground 

L  =  Litter 

Agcr  =  Agropyron  cristatum 

Saka  =  Salsola  kali 

Juos  =  Juniperus  osteosperma 


4. 


3 


TABLE  //  Tree,  shrub,  and  ground  cover  (percent)  on  windrowed  runoff  plot: 
at  the  Blanding  study  site,  September,  1971. 


Plot 


Transect    

Number  1/   Trees 


Percent  Cover 


Shrub 


Ground 


Windrow  #1    19  ft     0 


BG  58 

L  16 

Agcr  26 


33  ft 


55  ft 


7h   ft 


(Mean)  x 


Artr  3 


Artr  1 


BG  56 

L  lit 

Agcr  30 

BG  1+6 

L  29 

Agcr  25 

BG  65 

L  15 

Agcr  20 

BG  56 

L  19 

Agcr  25 


Windrow  #2    19  ft     0 


33  ft     0 


BG  52 

L  15 

Agcr  32 

Circium  spp   1 

BG  66 

L  15 

Agcr  18 

Sphaeralcea  spp  1 


55  ft     0 


7U  ft     0 


(Mean)  x 


BG  71 

L  k 

Agcr  25 

BG  68 

L  11 

Agcr  18 

Sphaeralcea  spp  3 

BG  61+ 

L  12 

Agcr  23 

Sphaeralcea  spp  1 


Table  //,  (contd) 


44 


Plot 


Transect    

Number  1/   Trees 


Percent  Cover 


Shrub 


Ground 


Windrow  #3    19  ft     0 


33  ft     0 


Artr  2 


55  ft     0 


7h   ft     0 


(Mean)  x 


BG 

57 

L 

27 

Agcr 

16 

BG 

1+9 

L 

32 

Agcr 

19 

BG 

1+5 

L 

39 

Agcr 

16 

BG 

i+o 

L 

1+7 

Agcr 

13 

BG 

1*8 

L 

36 

Agcr 

16 

BG 

5k 

L 

25 

Agcr 

21 

BG 

62 

L 

21 

Agcr 

17 

BG 

58 

L 

27 

Agcr 

15 

BG 

1+1+ 

L 

36 

Agcr 

20 

BG 

5k 

L 

28 

Agcr 

18 

Windrow  ttk         19  ft     0 


33  ft     0 


55  ft     0 


Ik   ft     0 


(Mean)  x 


Windrow  #5    19  ft     0 


33  ft     0 


BG 

1+2 

L 

32 

Agcr 

26 

BG 

65 

L 

11+ 

Agcr 

21 

Table  //   (contd) 


# 


Plot 


Transect  

Number  1/  Trees 


Percent  Cover 


Shrub 


Ground 


55  ft 


lk   ft 


.Mean)  x 


BG 

76 

L 

5 

Agcr 

18 

Sphaeracl 

ea  spp 

1 

BG 

69 

L 

5 

Agcr 

26 

BG 

63 

L 

lU 

Agcr 

23 

1/  Line  transects  across  runoff  plots  at  indicated  distances  measured  from 
top  of  plot. 


2/    3G  =  Bare  Ground 

L  =  Litter 

Agcr  =  Agropyron  cristatum 

EP  =  Erosion  Pavement 

Saka  =  Salsola  Kali 


k6 


APPENDIX  fll 
Soils  data  for  each  0.11-acre  plot 
at  Milford  and  Blanding 
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Abstract 

Lichen  crusts  on  soils  within  the  Colorado  Plateau  were  studied 
to  determine  their  effect  on  infiltration  rates,  potential  sediment 
production,  permeability,  and  several  chemical  properties  of  the  soil. 
Six  different  crust  stages  were  identified.  Undisturbed  soil  cores  were 
used  to  determine  intrinsic  permeability  under  three  treatments  and 
disturbed  soil  samples  were  analyzed  for  pH,  percent  organic  matter,  soil 
texture,  Ca+Mg  content  and  total  conductivity. 

It  was  found  that  the  lichen  crust  had  little  effect  on  soil  chemical 
properties.  Analysis  of  undisturbed  soil  core  data  indicates  that  high 
lichen  cover  tends  to  decrease  intrinsic  permeability;  this  effect 
was  reinforced  when  cores  were  irrigated.  Data  obtained  with  the  Rocky 
Mountain  infi Urometer  indicated  that  sites  with  any  degree  of  lichen 
cover  had  significantly  higher  infiltration  rates  than,  chained  areas  (no 
lichen  cover).  Patterns  of  sediment  production  indicate  a  potential  for 
increased  sediment  once  the  lichen  crust  has  been  disturbed. 


Introduction 

Research  has  shown  that  soil  crusts  formed  by  algae  and  other 
microflora  influence  infiltration  and  soil  stability.   Fletcher  and 
Martin  (6)  found  that  crusts  of  molds  and  algae  increase  tensile 
strength  of  soil,  reduce  erosion,  and  increase  organic  matter.  Booth  (2) 
reported  that  algal  crusts  resist  drought,  prevent  wind  erosion,  and 
serve  to  break  the  force  of  falling  raindrops.  Cameron  and  Blank  (4)  point 
out  the  importance  of  soil  crusts  in  preventing  erosion  and  in  colonizing 
bare  areas.  Soil  microflora  can  have  effects  on  soil  aggregation 
according  to  Bond  (1 ). 

The  objective  of  this  study  was  to  evaluate  the  influence  of  a  soil 
crust  composed  of  crustose  lichens  on  infiltration  rates,  potential 
sediment  production,  and  several  additional  soil  properties. 

The  lichen  crust,  as  found  in  the  Colorado  Plateau,  is  almost 
continuous  in  its  undisturbed  condition.  Many  small  lichen  "pedestals" 
about  three  inches  across  characterize  the  crust;  Lichen  species  found 
in  the  crust  include  Collema  coccophorum,  Derma tocarpon  hepaticum,  Fulgen 
si  a  fulgens,  Lecidea  decipiens,  Peccania  kansasa.  Lichen  identifications 
were  made  by  C.  W.  Wetmore.  These  species  are  predominately  b3ack  in  color 
giving  the  crust  a  dark,  "bumpy"  appearance. 

Procedures 
The  study  area  is  located  forty-five  miles  west  of  Blanding  in  San 
Juan  County,  Utah,  in  the  vicinity  of  Natural  Bridges  National  Monument. 
Six  conditions  of  the  crust  were  sampled: 

A.  relict  areas;  crust  undisturbed 

B.  grazed  areas  within  the  pinyon  (Pinus  edul us) -juniper  (Juniperus 
osteosperma)  woodland  (protected  since  1967)  where  the  crust 
was  still  intact  (near  relict) 


C.  grazed  areas  within  the  pinyon-juniper  woodland  (protected  since 
1967)  where  the  crust  was  in  an  intermediate  stage  of  breakdown. 

D.  grazed  areas  within  the  pinyon-juniper  woodland  (protected  since 
1967)  where  the  crust  had  been  severely  disturbed  (animal  pathways 
and  waterways) . 

E.  crust  mechanically  disturbed  by  double  chaining  (a  large  anchor 
chain  is  pulled  behind  two  tractors,  first  in  one  direction  and  then 
in  the  opposite  direction)  pinyon-juniper  and  leaving  the  debris 

in  place.  The  chaining  treatment  was  applied  in  the  fall  of  1967. 

F.  crust  mechanically  disturbed  by  chaining  the  pinyon-juniper 
woodland  and  windrowing  the  debris.   (This  is  a  common  range 
improvement  technique  used  in  the  area.)  The  chaining  treatment 
was  applied  in  the  fall  of  1967. 

Soil  samples  from  three  depths  (0-1/2",  1/2-1",  1-2")  were  taken  from 
sites  representing  each  crust  condition.  Soil  samples  were  analyzed  for 
the  following  factors  at  each  depth;  texture  (3),  degree  of  aggregation 
(Bouyoucos  (3),  with  Calgon  omitted),  Ca+Mg  present,  conductivity,  pH  and 
organic  matter  (7). 

Thirty-six  two-inch  diameter  soil  cores  were  obtained  from  each  crust 
condition.  These  cores  were  used  to  determine  intrinsic  permeabilities  (7). 
Water  was  ponded  over  each  core  and  a  constant  head  was  maintained.  Total 
percolation  through  each  core  was  measured  at  5  minute  intervals  for  a 
period  of  forty-five  minutes.  Before  permeability  runs  were  made,  one- third 
of  the  cores  from  each  site  were  subjected  to  one  of  the  following  treatments 

1.  cores  air  dry  prior  to  the  run 

2.  cores  wet  to  saturation  daily  for  a  period  of  ten  days  prior  to 


the  run  (cores  air  dry  for  10  days  prior  to  wetting  treatment) 
3.  cores  wet  to  saturation  every   five  days  for  a  period  of  two 

months  prior  to  runs 
Infiltration  and  sediment  production  measurements  were  made  for  each 
crust  condition  using  a  Rocky  Mountain  infi Urometer  (7).  Seventy- two 
infiltrometer  runs,  each  28  minutes  in  length,  were  made  (12  replications 
per  crust  condition).  Plots  were  prewet  and  soils  allowed  to  drain  to  field 
capacity  by  waiting  a  minimum  of  2  hours.  Pooled  sediment  production  samples 
were  taken,  i.e.  turbidity  of  runoff  was  averaged  over  all  time  periods. 
Infiltrometer  runs  were  made  during  the  first  2  weeks  of  June,  1971. 

Results  and  Discussion 
Average  values, for  all  measurements  regarding  each  crust  condition  are 
given  in  tables  1  /\  2,  and  3-. 

Soil  Properties 

Soil  properties  were  not  strongly  influenced  by  the  soil  crust.  Organic 
matter  (overall  average,  .42%)  and  pH  (overall  average,  7.29)  analyses 
showed  no  significant  differences  among  crust  conditions  or  soil  depths. 
Significant  differences  were  found  in  analysis  of  particle  size  distribution, 
Ca+Mg  content,  total  conductivity,  and  water  stable  aggregates  less  than  2 
mm.  diameter.  These  differences  can  be  attributed  to  differences  caused  by 
mechanical  disturbance  during  chaining  activities. 
a.  Texture 

Averaged  over  all  depths,  sand  percentages  were  significantly  higher 
in  the  chained  with  windrowing  treatment  (see  table  1).  Wind,  as  a  sorting 
agent,  is  more  effective  in  this  area  and  as  a  result  has  probably  deleted 


the  silt  fraction  slightly.  The  remaining  crust  conditions  did  not  show  this 
effect  since  their  surfaces  were  protected  by  lichen  cover  and/or  pinyon-juniper 
woodland  of  debris.  Silt  percentages,  averaged  over  all  depths,  were  like- 
wise significantly  different;  sites  A  and  B  were  significantly  higher  in 
silt  than  the  remaining  sites.  Clay  percentages,  averaged  over  all  depths, 
were  significantly  lower  in  sites  A  and  B  than  the  other  sites. 

The  clay  fraction  was  influenced  somewhat  by  soil  depth.  Averaged  over 
all  crust  conditions,  the  two  inch  depth  showed  significantly  greater  clay 
content  than  the  other  two  depths.  This  may  indicate  some  mobilization  and 
translocation  of  clays  to  lower  depths, 
b.  Physical  and  Chemical  Properties: 

Averaged  over  all  depths,  Ca+Mg  content  was  found  to  be  significantly 
higher  in  the  chain  with  windrowing  treatment  area  and  in  the  relict  stand 
(Site  A).  Total  conductivity  (averaged  over  all  depths)  was  also  significantly 
higher  in  the  chain  and  windrowed  area  than  in  the  remaining  sites. 

The  windrowing  treatment  involves  disturbance  and  mixing  of  the  soil 
to  a  depth  of  from  4  to  6  inches.  The  ion-rich  caliche  layer  which  is 
present  at  a  shallow  depth  in  this  geographic  region  was  disturbed  and  mixed 
with  the  surface  soil  during  mechanical  treatment;  this  probably  accounts 
for  the  higher  Ca+Mg  content  and  higher  total  ion  content  found  in  this 
crust  condition. 

Averaged  over  all  crust  conditions,  soil  depths  1.5  inches  and  2  inches 
were  significantly  higher  in  Ca+Mg  and  in  total  conductivity  than  was  depth 
.5  inch.  This  indicates  there  has  been  some  leaching  of  ions  since 
treatments  were  applied. 

Organic  matter  was  not  significantly  different  among  sites  or  depths. 


c.   Water  Stable  Aggregates: 

Water  stable  aggregates  less  than  2  mm  diameter  were  found  to  be 
significantly  higher  in  the  chain  with  windrowing  treatment  (averaged  over 
all  depths).  Averaged  over  all  crust  conditions,  depth  2   inches  also 
showed  a  significant  increase  in  stable  aggregates.  This  can  be  explained 
by  the  above-mentioned,  increase  in  total  ions  and  Ca+Mg  content  due  to 
mechanical  mixing  at  site  F  and  probable  leaching  phenomenon  at  depth  2 
inches. 

Intrinsic  Permeability 

Crust  condition  and  treatment  influenced  intrinsic  permeability.  High 
lichen  cover  significantly  impeded  permeability;  in  general,  the  greater 
the  lichen  cover,  the  greater  the  resistence  to  percolation.  Intrinsic  perme- 
ability did  not  vary  significantly  between  time  periods.  Significant 
differences  were  found  among  sites  with  appreciable  lichen  cover  (Sites  A,  B, 
C)  and  those  which  were  nearly  bare  of  lichen  cover  (see  table  2.).  Highest 
percolation  rates  occurred  within  cores  from  sites  with  yery   low  lichen 
cover;  lowest  values  occurred  in  cores  with  greatest  lichen  cover. 

The  two  irrigation  treatments  (see  Procedures)  significantly  decreased 
intrinsic  permeability  of  cores  as  compared  with  cores  not  irrigated. 
Addition  of  water  simulates  conditions  of  the  lichen  crust  during  favorable 
moisture  periods.  Irrigation  reinforces  the  tendency  of  high  lichen  cover 
to  impede  percolation. 

Infiltration  and  Sediment  Production 

Infil trometer  runs  indicated  that  high  lichen  cover  significantly 
increases  infiltration  rates.  Sites  with  any  degree  of  lichen  cover  (sites 
A,  B,  C)  showed  significantly  higher  infiltration  rates  as  compared  with 


Table  2.     Average  infiltration  rates   (cm/hr)  during  various  time  periods 
(minutes)  for  each  cryptogamic  crust  condition. 


Time  Period 

(Minutes)-/ 

Site 

3-8 

8-13 

13-18 

18-23 

23-28 

A 
(relict) 

9.3  ac 

8.1   a 

7.2  ac 

6.4  ac 

6.3  ac 

B 
(near  relict) 

10.5  a 

8.1   a 

6.8  ac 

5.4  be 

5.2  ac 

C 

(intermediate 
condition) 

8.8  be 

8.2  a 

7.6  a 

7.6  a 

6.7  a 

D 
(waterways) 

10.4  ac 

7.5  ac 

7.5  a 

6.2  ac 

6.0  ac 

E 

(chained-with- 
debris-in-place) 

8.2  bd 

6.3  bd 

5.8  b 

5.3  be 

5.0  be 

F 

(chained-with- 
windrowing) 

8.3  bd 

5.9  bd 

4.7  b 

4.6  bd 

3.8  bd 

—  Any  two  means  in  the  same  column  with  the  same  combination  of  letters 
are  not  significantly  different  at  the  .05  level  of  probability. 
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Doth  of  the  chained  areas  where  lichen  cover  has  largely  been  destroyed. 
There  were  no  significant  differences  among  unchained  sites  (A,  B,  C,  D) , 
but  all  were  significantly  different  from  the  chained  sites  (E,  F).  Lichen 
pedestals  help  prevent  surface  sealing  of  soil  pores  and  they  also  create 
some  detention  and  retention  storage  on  the  plots,  thus  increasing 
infiltration. 

Site  D  (pathways  and  waterways)  appears  to  be  a  special  case  in  regard 
to  infiltration  rate.  It  is  an  apparent  exception  to  the  generalization 
that  presence  of  any  degree  of  lichen  cover  enhances  infiltration.  There 
are  several  plausible  reasons  for  this  inconsistency.  Infiltration  measurements 
were  taken  in  the  spring  of  1971  when  the  frost  heaving  effects  of  the 
past  winter  were  still  yery   much  in  evidence.  In  this  arid  region,  soil 
moisture  accumulates  during  the  winter  months  and  resultant  frost  heaving 
causes  pronounced  "fluffing-up"  of  the  soil  and  thus  increased  friability 
and  decreased  bulk  density.  Past  experience  has  shown  that  increased 
infiltration  rates  are  to  be  expected  in  spring  and  early  summer  due  to 
these  effects.  With  this  in  mind,  one  could  speculate  that  once  summer 
thunderstorms  cause  flow  in  the  waterways,  the  apparently  high  spring 
infiltration  rate  would  be  reduced. 

Infiltration  rates  among  the  different  crust  conditions  remained  in 
the  same  positions  relative  to  one  another  throughout  the  infiltrometer 
runs  (e.g.,  a  site  with  higher  infiltration  rate  comparated  to  another  at 
the  beginning  of  the  infiltrometer  run  tended  also  to  have  a  higher  rate 
at  the  end  of  the  run;  see  fig.  2  and  Table  4).  Effect  of  yery   high  lichen 
cover  (sites  A  and  B)  on  infiltration  rates  was  greatest  early  in  the 
infiltrometer  runs  (see  fig.  2). 

Averaged  over  all  sites,  successive  time  periods  all  differed  significantly; 


the  first  five-minute  period  showed  highest  infiltration  rates  and  the  rates 
decreased  with  each  successive  period. 

Differences  in  sediment  production  among  the  crust  conditions  were  not 
significant  although  a  general  pattern  was  indicated.  The  trend  was  toward 
increased  sediment  once  the  crust  had  been  in  any  way  disturbed.  Degree 
of  disturbance  seemed  to  matter  little  in  regard  to  sediment  production 
(see  fig.  3). 

Conclusions 

1.  The  lichen  crust  has  little  direct  effect  on  the  soil  chemical  and 
physical  properties  considered. 

2.  High  lichen  cover  (sites  A,  B,  C)  caused  a  decrease  in  intrinsic 
permeability;  this  effect  was  reinforced  when  the  cores  were  irrigated. 

3.  Appreciable  lichen  cover  promotes  increased  infiltration  as  compared 
with  soil  surface  conditions  that  result  from  chaining. 

4.  Patterns  of  sediment  production  data  indicate  an  increase  in 
sediment  production  once  the  crust  has  been  disturbed. 
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Footnotes 
]_/     This  study  was  conducted  in  cooperation  with  the  Bureau  of  Land 
Management  under  contract  1*+- H- 0008- 283 7.   Journal  Paper  1241, 
Utah  Agricultural  Experiment  Station,  Utah  State  University,  Logan, 
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High) ight 
The  microclimatic  variables  net  radiation,  solar  radiation,  and 
wind  were  studied  during  parts  of  1968  and  1969  on  a  pinyon-jun iper 
site  in  southwestern  Utah.  Treatments  studied  included  cha in i ng-wi th- 
debris  windrowed,  cha in ing-wi th-debr is- in-place,  and  undisturbed  woodland. 
Net  radiation  was  greatest  over  the  natural  pi nyon-jun iper  stand;   Net 
radiation  on  the  cha in-wi th-windrow  treatment  and  chain-with- 
debr is- in-place  treatment  averaged  71  and  91  percent,  respectively,  of 
that  measured  on  control  plots.  Albedo  values  averaged  13  and  12  percent 
for  the  2  years  on  control  plots,  21  and  19  percent  on  the  chain-with- 
windrow  treatment,  and  13  and  ]k   percent  on  the  cha in-wi th-debr is- in-place 
treatment.   Roughly  3  miles  of  wind  (as  measured  at  approximately 
mid-canopy  height)  occurred  on  the  chained  treatments  for  every  1 
mile  measured  on  the  control  plots. 


The  pinyon-junipcr  (Pinus  spp.- Jun iperous  spp.)  type  extends  over 
approximately  60  million  acres  in  the  western  United  States  (Dortignac,  I960) . 
Of  that  amount,  several  hundred  thousand  acres  have  been  altered  through 
various  tree  removal  techniques  devised  by  the  Bureau  of  Land  Management 
and  the  Forest  Service.  The  trees  have  been  removed  by  chaining,  burning, 
bulldozing,  or  spraying  and  either  a  grass  species  was  artifically  seeded 
on  the  site  or  the  native  species  (shrubs  and  grasses)  were  allowed  to 
naturally  seed  and  eventually  occupy  the  site. 

Little  is  known  about  natural  or  undisturbed  pinyon- juni per  ecosystems, 
let  alone  those  systems  in  which  some  form  of  vegetation  manipulation  has 
taken  place.  This  study  was  designed  to  provide  some  preliminary  information 
regarding  microclimatic  changes  that  may  occur  when  pinyon-jun iper  sites 
are  significantly  altered  by  tree  removal.  The  study  site  is  located 
about  *+0  miles  southwest  of  Milford  in  southwestern  Utah. 

Methods  and  Site  Description 
The  study  site  was  chained  (using  a  large  anchor  chain  between  two 
tractors)  in  the  fall  of  1967.  Three  chaining  treatments  were  used;   chaining 
with  all  debris  windrowed  (crested  wheatgrass  (Agropyron  cr  istatum)  drill 
seeded  at  8  lbs. /acre);  double  chaining  with  debris- in-place  (crested 
wheatgrass  broadcast  seeded  at  8  lbs/acre);  and  no  chaining.  Approximately 
30  to  kO   acres  were  involved  in  each  treatment.  All  treatments  were  fenced 
to  exclude  livestock. 

Vegetation  on  the  control  plots  was  characterized  by  a  20  to  kS   percent 
canopy  coverage  of  pinyon  and  juniper  trees.  Shrub  cover  was  variable 
and  consisted  of  big  sagebrush  (Artemisia  trldentata),  black  sagebrush 
(Artemis  ia  nova) ,  and  small  rabbi tbrush  (Chrysothamnus  spp.).   Scattered 
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plants  of  Er  iogonum  spp. ,  Aster  spp. ,  and  Salsola  ka 1 i  were  also  present. 
Vegetation  on  the  cha i n-wi th-wi ndrowi ng  treatment  was  limited  during  the 
duration  of  this  study;  crested  wheatgrass  seedlings  were  the  dominant 
species  present.  Seedling  canopy  coverage  was  less  than  10  percent 
during  1968  and  from  10  to  20  percent  during  1969.   0n  the  chain-with- 
debr is- in-place  plot,  the  crested  wheatgrass  seedling  canopy  coverage  was 
less  than  10  percent  during  both  1968  and  1969,  ar,d  shrub  canopy 
coverage  (all  species)  was  15  percent  or  less  during  both  years. 

Soils  are  silty  loam  in  texture  and  were  derived  from  basaltic 
parent  mater  ials . 

During  the  summer  of  1 968  climatic  stations  were  established  in  each 
of  the  three  treatments.   Net  radiation  was  measured  at  20  feet  above  the 
ground  surface  in  the  undisturbed  and  cha ined-wi th-debr is- in-place 
treatments  and  at  8  feet  above  the  ground  surface  in  the  cha i ned-wi th- 
windrowing  treatment.   For  radiometers  mounted  20  feet  about  the  soil 
surface,  nine-tenths  of  the  field  of  view  is  within  a  diameter  of  115  feet 
(Reifsnyder  and  Lull,  1 965) .  A  Fritschen  (1963)  type  miniature  net 
radiometer,  as  modified  by  Campbell,  Ashcroft  and  Taylor  (196*+).  was 
used  to  make  the  net  radiation  measurements.   Measurements  were  recorded 
on  a  Rustrak  Model  88  12  v  dc  MA  recorder  with  an  unregulated  motor.   The 
radiometers  (and  all  other  instruments)  were  not  moved  about,  but  remained 
at  or  very  near  the  same  position  for  the  duration  of  the  study  (Figures  1, 
2,  3).   The  net  radiation  measured  is  the  difference  between  the  downward 
flux,  consisting  of  direct  and  diffuse  solar  radiation  and  atmospheric 
thermal  radiation,  and  the  upward  flux,  consisting  of  reflected  solar 
and  thermal  radiation  and  emitted  thermal  radiation. 

Incoming  and  reflected  solar  radiation  was  measured  with  a  Sol-A-Meter 


mm 


(silicon  cell)  from  the  Yellot  Solar  Energy  Laboratory.   Dirmhirm  (1968) 
has  discussed  the  use  of  this  instrument  in  radiation  studies.  A  single 
Sol-A-Meter  measured  incoming  solar  radiation  in  the  open  at  a  height 
of  approximately  6  feet.  All  measures  of  reflected  shortwave  radiation 
were  made  at  the  same  height  above  the  soil  surface  as  the  net  radiation 
measures.  Measurements  were  recorded  on  a  Rustrak  Model  88  12v  dc  MV 
recorder  with  an  unregulated  motor. 

Wind  was  measured  at  a  height  of  11  feet  (roughly  mid-canopy 
height)  with  Science  Associates  Model  992*+A  totalizing  anemometers. 

Results  and  Discussion 
net  radiation 

Net  radiation  values  for  each  chaining  treatment  are  shown  in 
Figure  k   for  the  period  8-24-68  to  9-22-68.   Data  for  the  control 
treatment  were  available  only  for  the  period  9-7-68  to  9-22-68. 

Highest  net  radiation  values  were  measured  over  the  control  treatment, 

followed  by  the  cha in-wi th-debr is- in-place  treatment.  The  chain-with— 

windrow  treatment  gave  the  lowest  net  radiation  values.   Net  radiation 

in  the  cha in-wi th-wind row  treatment  averaged  71  percent  (range  was  kS 

to  85  percent)  of  that  measured  on  the  control  plot.   Net  radiation 

on  the  cha in-wi th-debr is- in-place  treatment  averaged  9'  percent  (range 

was  80  to  105  percent)  of  that  measured  on  the  control  plot. 

The  above  results  are  somewhat  expected.   The  net  radiation  (R  )  of 

n 

a  surface  is  the  algebraic  sum  of  the  downward  flux  of  solar  radiation 
(Rs^)  from  sun  and  sky,  the  downward  infrared  or  thermal  radiation  flux 
(R j J)  from  the  atmosphere,  and  the  upward  flux  of  reflected  solar  radiation 
(Rsf) ,  and  the  upward  infrared  radiation  flux  (R\f)    from  the  surface. 


Thus: 

R   =  R  |  +  R.\  -    R  f  -  R.t  (1) 

n    s     i     s      i 

The  downward  fluxes  R  I     and  R . i     are  controlled  by  atmospheric 

si 

conditions  and  are  usually  rather  constant  except  under  partly  cloudy 

skies.   The  magnitude  of  R  f   is  dependent  on  the  reflective  properties 

of  the  surface,  the  ratio  of  outgoing  (reflected)  solar  radiation  to 

that  incoming  being  refered  to  as  the  albedo  (a).   Thus: 

R  f  =  aR  i  (2) 

s '      s  v  ' 

The  flux  of  R.f   is  dependent  primarily  on  radiation  emitted  by  the 
surface,  but  the  reflected  part  of  the  downward  infrared  radiation  is 
also  included.   Using  the  absorptivity  or  emissivity  e  ,  the  Stefan- 
Boltzmann  constant  <s  ,    and  the  fact  that  surface  temperature  (T  )  governs 
radiation  emitted  by  the  surface,  then 

R.t  =    drcrT^  +  (l-€  )  R.J,   .  (3) 

If  we  can  assume  that  €T  =  1  (usual  range  0.95  to  0.98  for  natural 

surfaces),  then  combining  and  rearranging  (1),  (2),  and  (3)  gives 

R   =  (1-a)  R  I    +  R.i  -  <rT  \  (k) 

n  s      1        o 

Examination  of  (k)   will  show  that  net  radiation  is  dependent  on 
the  albedo  and  surface  temperature  if  downward  components  (R  |  »  R-v  ) 
are  somewhat  constant.  These  conditions  were  closely  approximated  in 
this  study.   !t  will  be  seen  later  that  albedo  of  the  cha in-wi th-windrow 
treatment  was  considerably  higher  than  that  measured  on  the  other  two 
treatments.   Presence  of  debris  scattered  about  on  the  debris- in-place 
treatment  would  tend  to  reduce  the  albedo,  as  would  the  presence  of 
standing  trees. 

Influence  of  surface  temperature  differences  on  net  radiation  in  this 
study  is  not  known  at  this  time,  but  Federer  (1968)  has  found  in  the 


Northeastern  states  that  variations  in  surface  temperature  and  albedo 
were  about  equally  important  in  affecting  net  radiation.   Baumgartner 
(1965)  and  Tromble  and  Simanton  (1969)  have  found  that  net  radiation 
values  over  open  areas  may  be  only  about  60%  of  that  measured  over  a 
forest  canopy.  Tajchman  (1971)  found  that  total  net  radiation  over  a 
spruce  forest  in  Germany  was  20%  greater  than  that  over  alfalfa  and  16% 
greater  than  that  over  potatoes. 

solar  radiation 

Values  for  incoming  and  reflected  shortwave  are  shown  in  Figure  5. 
Data  were  available  for  parts  of  both  I968  and  I969.   During  the  I968 
measurement  period  the  albedo  of  the  cha in-wi th-windrow  treatment 
averaged  21  percent  (range  of  16  to  26  percent),  or  not  quite  twice  the 
albedo  of  the  cha in-wi th-debr is- in-place  treatment  (average  13  percent, 
range  7  to  20  percent)  and  control  treatment  (average  13  percent,  range 
9  to  21  percent) . 

The  above  relationship  held  for  the  I969  measurements  also.  Albedo 
of  the  cha in-wi th-windrow  treatment  averaged  19  percent  (range  of  12  to 
25  percent)  as  compared  with  ]k   percent  for  the  cha in-wi th-debr is- in-place 
treatment  (range  of  9  to  23  percent)  and  12  percent  for  the  control 
treatment  (range  of  k   to  22  percent).  Some  slight  changes  are  to  be 
expected  as  vegetation  increased  somewhat  on  both  of  the  chaining  treatments. 

Figure  6  gives  some  idea  of  the  magnitude  of  change  in  reflected 
solar  radiation  which  may  be  expected  following  a  snow  storm.   No  data  are 
available  for  incoming  solar  radiation.   The  increased  reflectivity  in 
response  to  snow  is  greatest  in  the  cha in-wi th-windrow  treatment,  intermediate 
in  the  cha in-wi th-debr is- in-place  treatment,  and  least  in  the  control 
treatment . 


lV-rf»»  ,-,//..  Ir*   ,ir/-~~)  , 


' lg^ i  i  1 1  i  i  i  i i m M ■ 


J***  s>«7s*s«**sa;jS5 


Figure 


i 

0 

* 

* 

*0 

ZSO- 

0 

1, 

\ 

B 

\ 

\ 

\ 

M 

S; 

Jtf  iso 

A 

A) 

^ 

$\/°°- 

& 

% 

$> 

^  •»■ 

* 

/7j 

**m******SW*MMMM*l8*fcU*«*** 


Figure 


Results  of  both  net  radiation  and  albedo  measurements  agree  closely 
with  those  of  Hornbeck  (1970)  who  compared  the  radiant  energy  of  clearcut 
and  forested  sites  in  West  Virginia. 

wi  nd 

Average  miles  of  wind  per  day  at  mid-canopy  height  as  averaged 
over  approximately  2-week  periods  are  shown  for  each  treatment  in  Table  1. 
The  relationships  are  similar  for  both  I968  and  1969.   In  general,  for 
each  mile  of  wind  recorded  on  the  control  plots  at  mid-canopy  height, 
approximately  3  miles  can  be  expected  on  the  two  chaining  treatments. 
Slightly  more  wind  was  recorded  during  all  time  intervals  on  the 
chain-wi th-windrow  treatment  than  on  the  chain-with-debr is- in-place 
treatment . 

Cone  1  us  ions 

The  objective  of  this  study  was  to  provide  some  preliminary 
information  regarding  microclimatic  changes  that  may  occur  when  a 
pinyon-juniper  site  is  significantly  altered  by  removal  of  trees. 
This  paper  has  provided  some  insight  regarding  net  radiation  differences, 
albedo  differences,  and  wind  relationships. 

Net  radiation  measurements  were  greatest  over  the  natural  pinyon- 
juniper  woodland.   Net  radiation  on  the  cha in-wi th-wi ndrow  treatment 
and  cha in-wi th-debr is- in-place  treatment  averaged  71  and  91  percent, 
respectively,  of  that  measured  on  the  control  plot.   Part  of  the 
difference  is  due  to  a  higher  albedo  on  the  cha in-wi th-windrow  treatment. 

Average  albedo  values  varied  only  slightly  during  measurement  periods 
in  I968  and  I969.  Albedo  values  averaged  13  and  12  percent,  respectively, 
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Net  radiation  measurements  were  greatest  over  the  natural  pinyon- 
juniper  woodland.   Net  radiation  on  the  cha in-wi th-wi nd row  treatment 
and  chain-wi th-debr is- in-place  treatment  averaged  71  and  91  percent, 
respectively,  of  that  measured  on  the  control  plot.   Part  of  the 
difference  is  due  to  a  higher  albedo  on  the  cha in-wi th-windrow  treatment. 

Average  albedo  values  varied  only  slightly  during  measurement  periods 
in  I968  and  I969.  Albedo  values  averaged  13  and  12  percent,  respectively, 


for  the  2  years  on  control  plots,  21  and  19  percent  on  the  chain-with- 
windrow  treatment,  and  13  and  ]k   percent  on  the  cha in-wi th-debr is- in-place 
treatment.   Presence  of  .debris  on  the  ground  in  the  debris- in-place 
treatment  tended  to  bring  albedo  values  close  to  those  measured  on  control 
plots,  at  least  during  snow-free  periods.   Following  a  snow  storm,  the 
amount  of  reflected  solar  radiation  increased  most  on  the  chain-with- 
windrow  treatment  and  least  on  the  control  treatment. 

As  might  be  expected,  wind  (as  measured  at  approximately  mid-canopy 
height)  was  greatest  where  standing  trees  were  absent.   Roughly  3  miles 
of  wind  occurred  on  the  chained  treatments  for  every  1  mile  measured 
on  the  control  plots. 


Table  1.   Summary  of  wind  measurements  at  Milford  Study  Site. 


Time  Interval 

< 

Number 

of 
Days 

Average  Mi 
Control 

les  of  Wind 

Cha  in, 

Windrow 

per  Day  — 
Cha  in 
Debris- in-Place 

*7-  14-68  to  7-27-68 

13 

48 

144 

118 

7-27-58  to  8-9-68 

13 

32 

111 

97 

8-9-68  to  8-23-68 

14 

75 

210 

175 

8-23-68  to  9-6-68 

14 

42 

142 

137 

9-6-68  to  9-20-68 

14 

60 

181 

163 

9-20-68  to  10-3-68 

13 

38 

148 

142 

10-3-68  to  10-18-68 

15 

51 

160 

144 

i 
10-18-68  to  11-1-68 

14 

30 

119 

111 

.1 1-1-68  to  11-14-68 

13 

40 

149 

136 

11-14-68  to  11-30-68 

16 

34 

136 

121 

4-4-69  to  4-21-69 

17 

45 

— 

141 

,4-21-69  to  5-2-69 

11 

64 

195 

174 

5-2-69  to  5-17-69 

15 

41 

156 

144 

5-17-69  to  5-30-69 

13 

51 

166 

153 

5-30-69  to  6-12-69 

14 

43 

140 

117 

6-12-69  to  6-25-69 

13 

36 

131 

119 

6-25-69  to  7-11-69 

16 

61 

185 

161 

•*   At  11.5  feet  (approximately  mid-canopy  height)  above  the  ground  surface 
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Figure  Ti  ties 

Fig.  I.    Instrumentation  on  control  plot  in  the  pinyon-juniper  type. 
Fig.  2.    Instrumentation  on  cha in-wi th-debr is- in-place  treatment. 
Fig.  3.    Instrumentation  on  chain-wi th-windrow  treatment. 
Fig.  *t.   Net  radiation  as  measured  on  various  dates  over  each 

treatment . 
Fig.  5.    Incoming  and  reflected  solar  radiation. 
Fig.  6.   Reflected  solar  radiation  on  the  three  treatments  following 

two  periods  of  snow.  Squares  represent  the  windrow  treatment,  circles 

the  control  treatment,  and  dots  the  chain-wi th-debr is- in-place 

treatment. 


Bur.au  of  Land  Mam***,** 

Library 

Denver  Service  Center 


/ 

USDI -  BLM 

DATE 
LOANED 

QK 
495 
jg      Influence  of  Chaining 

G53       Watershed  Values  in  Ut 
1972 

Co 

00  ^ 

u 

CD 
O 
XI 
XI 

o 

m 

X 

n 

0 
O 
W 

50 
O 

PI 

50 

